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ABSTRACT 


For infrared sensors located at the nose of a missile or 
a projectile, an ace-olda problem occurs. A conflict between 
optical quality and aerodynamic drag exists. A new approach 
for solving this problem is being developed using gradient 
refractive index, GRIN, for lens design. The equations for 
lens design incorporating GRIN are formulated for the special 
case of spherically symmetrical GRIN. A computer code for 
designing a GRIN lens and for determining aberrations has 
been developed. A different program (LENS) calculates lens 
design and lens performance for the case of homogeneous 
refractive index. Results of this program provide a check 
on the more complex program which includes GRIN. The 
computed program calculates the spot diagram for both mer- 


idianal and skew rays. The spot diagram, which is the inter- 


section of rays with the focal plane, is plotted, and the Bs 


ms 
and centroid location of the spot diagram are calculated. 


Also a cumulative energy diagram is obtained from the progran. 
The computer code has been organized to calculate the items 
enumerated above for either a circular cone for the inside or 
outside surface of the pointed lens. For homogeneous 
refractive index, both cases are included in the program 
discussed in this report. However, for the GRIN case, only 
the complete program for a lens with a conical surface on 

the outside is reported here. The alternate case is being 


developed by a thesis student. 
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Features of the computer program include the following: 


use of many subroutines for clarity and for easy modification, 


and use of program structure for easy integration into an 
optimization code. 

Sample lens designs and lens performance are presented. 
The computer code provides the tools for lens design. Consider- 
able additional work is needed for selecting the best lens 
design. 
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1. INTRODUCTION 


‘ A long-term conflict exists between optics and aero- 
dynamics for missile nose design. Optics requires a blunt 
shaped nose similar to that used on the Sidewinder missile. 
Favorable aerodynamics for supersonic missiles requires a sharp 


mose which permits an attached shock wave. For rocket 


powered missiles or projectiles, the sensor and associated 
optics occupies the complete nose area. The penalty of a blunt 
nose is high drag. For ramjet propelled missiles or pro- 
jectiles, the penalty of a blunt nose is poor pressure recovery 
by the inlet; this statement assumes a nose inlet is used. 

In times past, the lens designer had available three design 
parameters; a satisfactory lens in the general shape of a cone 


has not been designed. A lens can be designed using homo- 


geneous refractive index which focuses rays; however, off-axis 
rays have very large aberrations. The thesis by Terrell [1] 
discusses the design of a conical lens. 

A fourth design parameter has become available which is 
spatially variable refractive index. The associated technology 
is known as GRIN; see references [2-13]. GRIN is short notation 
for Gradient Refractive Index. The thesis by Frazier [14] 
briefly introduced the use of GRIN in lens design; likewise 
the thesis of Terrell [1] has a brief discussion. The thesis 
of Carr [15] will have more extensive results. 


The development of materials which have variable refractive 
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index has been motivated by optical fiber technology. Optical 


fibers can propagate a single mode with low losses if an ap- 


propriate spatial variation of refractive index exists. Based 
on the information in reference (2] , the maximum variation 
of GRIN is limited to a few percent. The report contains 
twelve sections. In section 3 the problem of the homogeneous 
case is solved (program LENS) while the GRIN lens design is 
formulated in section 4 (program GRIN). The flow charts to 
the computer programs are given in section 5 and 6, and a 
user's guide to each of the programs is illustrated in section 
7. The programs are listed in sections 9 and 10. 
2. BRIEF OVERVIEW OF PROBLEM FORMULATION 

Define the axial distance as x and radial distance as 
r ina cylindrical coordinate system. The lens is to be 
axisymmatric so that angle is not a variable. 

In general the shape of the front surface, is a function 
of fr 

x = F,(r) (2.1) 
where subscript f£ denotes front surface. Likewise the rear 
surface is a function of r 
x = FL (r) (2.2) 

where subscript r indicates the rear surface. 

As a simplification, either the front or rear surface is 
a circular cone. Terrell [1] treated the case of a conical 
front surface; Carr [15] is investigating the case of a conical 
tear surface. This report develops a computer code which 


calculates either case. A single consistent computer program 


is used for both cases; this fact is important for application 
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of optimization codes. 
In general, the GRIN can be specified as 


n = F(x,r) (2.3) 


The function F involves both spatial variables. However, a 


less complex GRIN is used which is spherically symmetrical 


refractive index. In this case 


n= n(x ,r') (2.4) 


where r' is radial distance from the center of symmetry for 


the sphere. 


The location of the center of the spherical GRIN 
is on the lens axis at axial position Xo: 
One very important computational advantage of spherical 


GRIN is that the rays within the lens are confined to a single 


plane. In the case of skew rays, the plane of the rays is 


not a meridianal plane of the lens. 


The equation for a ray trajectory in a spherically 


symmetric GRIN geometry is 


dr' 
Cae 8 + =| 3 (2.5) 


r'“n° (r') -E 


where E is a constant for a ray. The value of E is obtained 


from 


E = ENy (rg) rgsin¥y (2.6) 


Define r' as a unit vector in the direction of a line extending 


from x, to point P in Figure 1. Define R as a unit vector of 


a ray which has been refracted by the surface at point P. The 


angle % is obtained from 
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Figure l. 


cos ¥) = reek (2.7) 


CENTER OF SYMMETRY OF GRIN 


Geometry for Ray Tracing in a Spherically Symmetrical 
GRIN medium. The ray vector, , and the radius 
vector, r' , are illustrated. 


In equation (2.6), the e€ is a sign function having a value of 


+ 1.0 or - 1.0 depending on the value of ¥ within the lens. 


+ 1.0 Y < 1/2 
€ -| (2.8) 


- 1.0 Y > 1/2 


Equation (2.5) can be integrated for certain functions n(r'). 


One function is 


n(r') = |A + Bir'/r)* 5 (2.9) 


where A and B are constants. rp) is a value associated with a 
specific lens design; a convenient value for Xo is that value 
for the ray defining the lens aperture. 
Combining equations (2.5) and (2.9) and completing the 
integration yields 
yore 
6, (r") = -k | sin or (sin |r 
(A -+4BE /%q ) (A°+4BE [Xo ) 
2.10) 
Equation (2.10) is an explicit relation for ®, as a function 
of r‘. res U, are the radii to the rear surface and to the 


front surface respectively. An explicit relation for r' as a 


function of ee can be derived: the result is 


73 |E| 
rr, f) (6,) = bs 
aa A+ A°+4BE Jo" * sin +26 +sin 
; to pa aa 2 
AY+4BE“/rp 
{ (2.11) 
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For details look at section 4.1. For certain phases of the 
calculation equation (2.11) is needed. 

Thinking ahead to numerical optimization, the function 
of equation (2.9) was chosen to avoid numerical integration 
of equation (2.5). Equation (2.5) must be evaluated several 
times for each ray. A lens design may involve 100 to 1000 
rays. Numerical optimization of the lens requires calculation 
of many different lens shapes. Hence, evaluation of equation 
(2.5) may occur several thousand times in the process of lens 
optimization. The GRIN-related parameters to be determined 
in optimization are A, B and Xo- 


3. EQUATIONS FOR THE HOMOGENEOUS CASE 


3.1 Meridian Plane 
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The homogeneous case is included here for several reasons. 
First, the homogeneous case is a simple example in which the 
intricacies of lens design are apparent in considerable clarity. 
Iteration is not required for design of a lens with a homo- 
geneous refractive index, whereas for GRIN, iteration is 
Necessary. Second, the homogeneous case provides a meaningful 
check of the GRIN program. To conduct the check, the value of 


B is zero in equation (2.9). Obviously the answers from GRIN 


program with B equal to zero must agree with.answers from the 
homogeneous case. 


The book by Kingslake [16] has been helpful in the program 


development. 
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3.1.1 The Equations 


: Refer to Figure 3 for definition of the symbols. 


n(J) + T,' (9) = 8(TI) 
sin(I,[J]) = N32 ° sin(I,"(J]}) 
The notation N39 is used for the ratio n;/n2- 


I,°(9) = (t/2-a) = (n(J) + I, (3)) 
sin(I,[J]) =n ° sin(I,'(J]) 
@+I1j,(J) +U = 1/2 


3.1.2 Solution 
The angles a, &8(J) are inputs and are known. 
From eq. (3.5): 
T, (9) =nf/2-a-0 
From eq. (3.4): 
I,'(J) = sin *{n, > sin(I,[J]) ] 
Substitute eq. (3.5) to eq. (3.3): 
n(J) + T5 (9) = T, (9) - I,‘ (J) + U 
But from eq. (3.2): 
I5'(J) = sin” *{ng. ° sin(Z5(J1)] 
From eq. (3.3.1): 
T4(3) = (1,(9) = 1,°(9)] = n(3) <0 
Substitute eq. (3.3.2) into eq. (3.2.1): 
T,'(3) = sin*In, * sin{I,(3) - 1,'(3) - n(3) + U}] 
Substitute eq. (3.2.2) into eq. (3.1): 


(3.1) 
(3.2) 


(3.3) 


(3.4) 


(3.5) 


(3.5.1) 


(3.4.1) 


(3.3.1) 


(3.2.1) 


(3.3.2) 


(3.2.2) 


n(J) + sin”*[ny, * sin{z, (3) = Ty'(3) - n(g) + UF] = E(3) 


(3.1.1) 


a ge 


Equation (3.1.1) is solved for n(J): 
223 


Use the formula: 

sin(A + B) = sinA cosB + cosA sinB 
And substitute: 

x = sin(n(J)) ; “7 _ 2 = cos(n(J)) 


1l1-x 
Equation (3.1.2) will turn into the form: 
1 
x = 
Aer 
Where: 


N53 ° cos [I, (J) - T,'() + U) = cos[B(J)]) 
A= 
No3 ° sin{I,(J) - I,'(3) + U) - sin(6(J)] 


Therefore: 


1 


n(J) = sin ~(x] 


Knowing n(J), I5' (J) can be solved from equation (3.1), and 


then, I,(J) can be solved from equation (3.2). 


. sin (I, (J) - I,' (J) = n(J) + U) = sin[{ AJ) - n(J)]j 


. RUE 6 


(3.1.2) 


(3.1.3) 


(3.1.4) 


(3.1.5) 


3.1.3 Coordinates of the Front and of the Rear Surfaces 


3.1.3.1 The First Ray 


All the magnitudes will be normalized 


with respect to the radius R. R is the aperture for the lens. 


Therefore: 
R = 1.0 


R 
x (1) © gente 


10 


(3.6) 


The thickness of the lens is defined: 


n(1) = fxg) = x,(1? + ty,() = y,( 1? (3.7) 


Therefore, the appropriate coordinates of the rear surface are: 


Xo (1) = x, (1) + T(1) + cos[n/2 -a - T,°(1)] 
(3.8) 
Yo(l) = y, (1) - T(1) + sin[1/2 - a - T,'())] 


3.1.3.2 The Other Rays: Coordinates of the Rear 


Surface 
The new point of the rear surface 
x5(J),y5(J) will be the intersection between the surface and 
the ray. 
; (a) The equation of the surface: 
Y2(J) = m(x5(J) - BB] (3.9) 
m= tan[(t/2 - n(J - 1)) = cot{n(J - 1)] (3.9.1) 
BB' = BF - B'F 
B'P = yo(J - 1) * [cot{"/2 - n(J - 1)] + cot{B(s - 1))) 


ee B'F = yo(5 - 1): {tan[n(J - 1)] + cot(B(s - 1) J} (3.9.2) 


ee ¥2(J) = cot[n(J = 1) }{x5(J) - [BF - B'F]} (3.9. 3) 
(b) The equation of the ray: 
Y2(J) = ~tan(B(J)]x5(J) + BF + tan [8(J)] (3.10) 
eeyg(J) = tan[B(J)] + {-x,(J) + BF} (3.10.1) 


From equations (3.9.3) and (3.10.1), the coordinates of the 


next point at the rear surface are found; those are: 


x.(3) = BR{tan[8(o)] + sot (aig -_1)]} - cot{[n(J-1)] > B'F 
2 cotin(j = + tan 


¥2(3) = tan[8(J) }{-x,(J) + BF} (3.11) 


eo Seamed 


li 


3.1.3.3 Coordinates of the Front Surface 


(a) Equation of the rear surface: 


¥,(9) = (tana) ° x, (J) (3.12) 
(b) Equation of the ray: 
¥,(J) = m[x (J) - x5(9)] + yo(d) (3.13) 
where: m = -tan[{t/2 - a - I,'(J)] (3.13.1) 
or: m= -cot[a + I,‘ (3) J (3.13.2) : 


From equations (3.12) and (3.13) the intercept point is 


calculated: 
Y2 (J) + Xo (J) * cot{a + T,'(9)] 


x, (9) = 
tana + cot[a + 1,'(J)) 
(3.14) 


y, 9) = (tana) - x1 (9) 
3.1.3.4 The Thickness at the New Point 
n(a) = / (x,(3) - x,(3)1? + fy, 0) - y5(9)1? (3.15) 
ea 1 Yi Y2 : 


3.2 Skew Rays in the Homogeneous Case 


3.2.1 Coordinate Transformation 


The x,y coordinates will be shifted by the 


angle a5- 
' = e - e e ° 
x x cos (a.) y sin (a) (3.16.1) 
y' =x: sin (c.,,) +y° cos (a,,) (3.16.2) 
z's Zz (3.16.3) 
3.2.2 Calculate the Intercept Point of the Front 
‘ Surface 


Assuming that the (0,0) point is at point A 


(Refer to figure 2.), the following equations can be derived: 


(a) Ray: Wie ay? + 22 (3.17.1) 
Yo! 2 2 
°, = - e t 
ee Ray costa, CS Be) tanto) * 29 
(3.17.2) 
Where Reay is the radius of the ray and the values for Yo's 
Zo" are known. 
(b) Cone: Reoné = tan(a) ° [% + AB] (3.18) 
Solving equations (3.17.2), (3.18) for x results: 
CDS” + A.D, * CZD 
XO = = ,DS 5 22 (3.19.1) 
2°2 2°2 
Y,' 
ZO = Zo! (3.19. 3) 
Where: 
A.D, = (tana)? = (tana_)? (3.20.1) 
2 2 2) e e 
Yo' 
CDS = cosa, tance, + (AB) [A5D,] (3.20.2) 
2 ' 
Yo' 2 2 2 Yo 
CZD = (seats + %)'° = (AB) “[A,D.]* - 2 (AB) soaa tana, 
(3.20.3) 
3.2.3 Direction Cosines of the Normal 
¢ ay? + z = (tana)?(x+ AB)? = 0 (3.22) 


éf 


st. -2(tana) 27(x + AB) (3.22.1) 


Figure 4. Diagram showing coordinate transformation for 
off-axis rays. 


is = ; SE = 22 (3.21.2) 
[Af] = 2¢(tana)4(x + aB)? + y? + 227% (3.21.4) 


Substituting equation (3.21) can tum equation (3.21.4) into 


the form: 
= 2tana 
lane] <aas (x + AB) (3.21.5) 


The direction cosines of the normal are therefore: 


ry aw SSELEX) oping (3.22) 
lat] 
(8£/Sy) 2y 
ga ELOY) ee (3.23) 
"Y lat] 2tana/cosa (X + AB) 
Define: o = tan~2(2z/y) (3.23.2) 


y? + 2 » (tana) 2 (x + AB) * 
y7(1 + tan2c) = tan2a(% + AB)? 


y = (coso) (tana) (x + AB) (3.23.2) 


Therefore: 


ny = (cosa) (cosa) (3.23.3) 
Similiarly: 
ny = ae = (sinc) (cosa) (3.24) 


3.2.4 Refraction at the Front Surface 


‘ The skew ray R will be defined: 
R= B+ () +R * (&) +R + (ep) (3.28) 


15 
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Figure 5. Cross section of lens in a plane normal to the lens 
axis. Point P is the point where an off-axis 
ray intersects the front surface of the lens. 
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Where Ry: Rys Ry are the direction cosines of the ray, 


and they have the value: 


Ry = cosa, i R= ~sing, + R, = 0 (3.26) 
The surface normal is defined by the following expression: 

> aA aA “a 

N=ny° (ey) + ny. (ey) + ny (e,) (3.27) 
Where the direction cosines of the normal surface were 


defined previously (equations 3.22, 3.23.3, 3.24). 


The refracted angle at the front surface is defined: 


> 
@ = cos? RN (3.28) 
IR]. Isl 
Where: |R| = JR,” + < + oe =1 
(3.29) 
i ee 2 2. 
There fore: 
= -1 
8 = cos {Ry ~My +R. ny tR, - nz} (3.30) 
¢= 7 = 8 
“iy (3.31) 
¢' = sin [N15 - Sind] 
3.2.5 Refraction Equations 
Referring to figure 6 the refraction at the 
rear surface with regard to the front surface will be as 
follows [16, p. 147]: 
nok! = n,R + [n, cos (¢") - n,cos(¢) JN 
The direction cosines of the normal surface were defined in 
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equations (3.22), (3.23.3), (3.24) and the direction 
cosines of the ray at the front surface were defined in 
equation (3.26). The direction cosines of the skew ray at 


the rear surface, will therefore be: 


K' = N15 R, + {cos(¢') = Ni cos(¢)] . ny 


L' = N15 Ry + [cos(¢') - Nj> cos(¢)] . ny (3.32) 


= 
u 


N12 R, + [cos (¢') - Ni> cos(¢)] . ne 


3.2.6 Intercept at the Rear Surface 
The equation of the cone of the rear 


surface is: 


y? +27 = [x -b)? (3.33) 


where a is the slope and is defined from figure 7. 


R R - R 


J JI+i1 rey 
a= tan(—&.) = — = _ : (3.33.1) 
i i aie Se ao 


b is the intercept at the X - axis and is calculated from 


equation (3.33.1): 


R e x ” R ° Xx 
bex, = Sti 2 jt (3.33.2) 
J+1 J 
(XI, YI, ZI) is the point of intercept at the rear surface. 
The equation of the rear surface is therefore: 
2 2 
YI" + ZI” - (aXI - b) = 0 (3.33.3) 


Where a,b were cefined in equations (3.33.1), (3.33.2). 


Figure 6. Diagram used to calculate the direction cosines 
of skew rays at the rear surface. This diagram is 
based on Kingslake [16]. 
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Figure 7. Geometry of a truncated cone which represents a 
small portion of the rear surface of the lens. 
The cone angle is er 
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The direction cosines of the ray at the rear surface were 


defined previously (3.32). They are: 


+. XI = XO , -, YI - YO. ~ 21-20 
K =z =p) ; L =z as » rece ? M res 5 REE (3.34) 
(xO, YO, ZO) is the point of intercept at the front surface 
(3.19.1, 3.19.2, 3.19.3). 
———— 
2 2 2 
D = /(XI = xO)“ + (YI - YO)“ + (ZI - 20) (3.35) 
Solving for D from equations (3.33.3) and (3.34) one gets: 
Page ° vp 2 - a'c' 
D ar + er (3.36) 
Where: 
at = (L')2 + (mt)? = a@(xry? (3.36.1) 
b' = (L')¥O + (M')ZO = (K')a*(XO = b/a) (3.36.2) 
c' = yO" + 20% - a*x0(xO = 2b/a) - b* (3.36.3) 
Knowing D, the intercept point at the rear surface can be 
verified using equation (3.34): 
XI = K'D + XO (3.37.1) 
YI = L'D + YO (3.37.2) 
ZI = M'D + ZO (3.37.3) 


3.2.7 Equation of the Normal to the Rear Surface 


The equation of the cone of the rear surface 
4 was defined previously (eq. 3.33): 


£(X,Y,2,) = 7(ax -b)? +y% +27 =0 (3.33) 
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ore ag ies a 


a,b are defined in equations (3.33.1), (3.33.2). 


Bee eae Leroy 
Similar to what was done at the front surface: 
es = ~2a(ax - b) (3.38.2) 
Sf. 2 , £22 


Pel= chr? + (8h) 24 Sh 
= 2[a*(ax - b)? + e + 2717 
= 2(ax - b)[(1 + a*) 4] (3.38.4) 


The direction cosines of the normal to the surface are as 


follows: 

iS-£ 8X -a ( 

ne ee MELE). 3.39.1) 

4 fe (1+ a’)? 

ny = pay (3.39.2) 
(y" + 2°) “(1 + a”) 

- z 

Tega eRe se ewe +39. 
By? + 2°) F1 + at) soeanve! 


3.2.8 Refraction at the Rear Surface 
Similar to what was done at the front surface: 
Ri = K'(e,) + L'(e,) + M' (ey) (3.40.1) 
N = ny (ey) + ny (ey) + ny (ey) (3.40.2) 
Where the direction cosines of the refracted ray (R') and 


the normal (N’) were defined previously (equations 3.32, 3.39). 


The refracted angle at the rear surface: 


> 
- N 
| 


NT! 


Both the ray and the normal are normalized; therefore: 
oy = K'ny + L'ny + M'n, (3.41.1) 
- =l1 
oy' = sin [n53 ; sin¢,] (3.41.2) 
Direction Cosines of Rear External Refracted 


Ray 
As done previously (equation 3.32): 


K" = No3K' - {cos ($,") - ny 3008 ($7) ] » ny 
L" = No3L' - (cos ($,") - n5,008 ($7) ] » Ny (3.42) 


zZ"* = No32' - [cos (¢,') ~- Ny 3Cc08 (7) ] - Ney 


3.2.10 Intersection with the Image Plane 
Xt = BF = D3K" + XI 


A 
Yim = DL + YI 


= nw 
Zim D3M + ZI 


Solving equation (3.43), the intercept with the image plane 
can be calculated: 
D3 = (BF = XI) /K" 
Y= (BE 5X) L" + yI (3.43.2) 
7 (ys XI )M" + 2z 
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4. FORMULATION OF EQUATIONS FOR GRIN LENS DESIGN 


4.1 Solving the Integral for a. 


As mentioned in section 2: 


x 
a 
5 r 


@ = 6, +€E ee 
= 0 r{r“n(r)“ - E*} 
ad 


te, TF, are the radii to the front surface and to the rear 


surface, respectively. 


Eze. No (F¢) - Xp - sinv, 


n = A+ B(E) 


Substituting: r 
x 

0 dv. _ - 

ve(s) oo: gps -2zq°r 


Equation (4.1) turns into the form: 


Vr 
eee ee dv 
e 0 284 2 EB /% 
Vv jav +b-V = 
£ Lo 
Therefore: (X9/r_) 


2. 
_7 (2(E/rt,)“V - A 
8, 7-2 | sin > pe 
(A“” + 4BE [Xp ) 2 
(9/T_¢) 


or: 


A” + 4BE /Xo ) A“* + 4BE /tq ) 
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(4.1) 


(4.2) 


(4.3) 


(4.1.1) 


(4.1.2) 


2 2 
2E - A 2E"/ - A 
a ea rie tgs eka 
: ( 


The appropriate expression for the radius will therefore be: 


Vv 2/B! 


ae re er 


/ 


28° /r oe 2) . 


Ja? + 4BE*/r,° 


(4.4) 


» 


Va + Vat + 4BE*/x,7 


- sin{#20, + sin” 


If the radius to the rear surface (r,) is calculated from a 


known radius to the front surface (r¢) the minus sign is used 
in front of the angle @,. This situation occurs with skew 
ray calculation. The plus sign is used when ie is calculated 
from a known rh: This situation occurs for lens design. 


4.2 


Meridian Plane 


4.2.1 Solution for the First Ray 
(a) Define Point X, (1) ,¥, (2) 


As in the homogeneous case: 


R= 1.0 


¥,(1) = R 


R 
%,(1) = fang 


Define Angle to the Front Radius 


Re ea. 1 SR Tee 
sin (6) (1) + 6, (1)) 


Refering to figure 8, 85 is the angle to the rear radius, 8 


is the angle between the rear and the front radii, and Xa 


is the location of the center of symmetry. 


{c) Choose 8. from Initial Thickness 


aL) = TCL) /r,(1) 
(4.8) 
Qo(L) = (8p(1) + O(L)) = 0602) 


(d) Solve for Ve 


a + Tj (2) + us 1/2 
(4.9) 


T, (2) =m w1T/f2-a-u 


Knowing r,(1), the initial index of refraction (n5(r_)) 


can be calculated from eq. (4.3). Therefore: 


sin(I,(1)) = nylrg)/my - sin(T,"(L)) 
1,'(1) = sin”*{nj/ng(r,) . sin{l, (1) 1} (4.10) 


Refering to figure 8, ve can be calculated. 
8,00) + 8y(1) sat 6 
oS 8, (2) + 84 (2) -@ (4.11) 
6 + bell) - I,‘(1) = 1/2 
Well) = m/2 + T'(1) ~ 6 (4.12) 


(e) Solve for £ 


Ese. n (kg) ~¥e- sinv, (4.2) 
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JOW JON 


eo a ee 


(4.2.1) 


ys 1/2 +e = +1 
Where: e = tl] 


y > T/2>+>€ 2 -1 


(£) Solve for rer¥, 


Knowing the radius to the front surface 
(r¢) and the angle between the two radii (6), the radius to 
the rear surface (x, can be calculated using equation (4.4). 


Using equation (4.2), vy, can be calculated: 


sol E 
vy = sin fs Holz) >», ty (4.2.2) 


r 


(g) Solve for 8(1) 
The focal length, which is an input, is known: 
BF = F/R 


The coordinates at the rear surface are: 


¥, (1) = r,sin(@)(1)) 
(4.13) 
X5(1) = Xo + r,cos (8, (1)) 
Using equation (4.13), 8(1) can be calculated: 
Y,(1) 
(1) = can? | = | (4.14) 
2 


(h) Solve for the Refractive Angles at the 


Front Surface 


As in the homogeneous case: 
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ane ata ne se a 


a a a eT 


n(J) + I,‘ (J) = B(J) 


sin[I,(J)] = n3/no (rz) e sin(I5'(J)] 


v (9) + I, (J) + n(J) + Oy (od) = 7 


Solve for n(J) in a manner similar to the homogeneous case: 


-1] / 1 
n(J) = sin (sy | 
V1l+A 
Where: cosé' - no/no(r.) . cos8 (J) 


siné' - nj/no(r,) - Sing (J) 


Where: 
6 = - 94 (9) = Yo (9) 
Knowing n(J), T5" (9) can be solved from eq. (4.15) and then, 
I, (a) can be solved from equation (4.16). 
4.2.2 Solution for Next Rays 
4.2.2.1 Coordinates at the Rear Surface 
Solve for the subsequent rays the 

same way as in the homogeneous case (equations 3.10, 3.10.1) 


4.2.2.2 Calculate Radius and Angle ’.to ‘the 


Rear Surface 


/ 2 2 
ri (J) a, (X4 (J) - Xo) + Yo(9) 
-1 
89 (9) = tan [Yo (J) /(X, (d) - X_) ] 
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(4.15) 


(4.16) 


(4.17) 


(4.18) 


{ 


(4.18.1) 


(4.18.2) 


(4.19) 


(4.20) 


4.2.2.3 Solve Front Surface Parameters 
Refer to Figure 8 for a definition 
of geometrical variables. 
(a) Equations 
n(J) + I4'(J) = B(d) 
sin[I,(J)] = n3/no(r,) : sin[I,'(J)] 
v (9) + I, (J) + n(J) + 6,(J) = 7 
BS e:%< no (xr) a sin[p (J) ] 
E=e. ny (r¢) - il, - sin[ (J) ] 
O,(5) + 89(J) sa + 6 
6 + (do) - Ty") = 1/2 
re = f(r, 6,(9)) 
sin{I,(J)] = No (re) /ny - sin(I,‘(g)] 
a + I, (9) +us= 7/2 


(b) Solution by Iteration 


An explicit solution is not 
possible. An iterative solution is necessary. 
(b1) Assume @, (J). 
(b2) Soive for coordinates 
of the front surface: 
- Calculate 1, (v) from eq. (4.24) 
~ Calculate 6 from eq. (4.21) 


- Calculate front surface coordinates and radius: 
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| 
| ¥, W) = tan(8) (J) + 8, (9)) . (X, (J) - Xo) (4.25) 
¥, (9) = tana. X, (9) (4. 26) 
The result is: 
X, . tan[6,(J) + 6 (J)] 
a Cc 0 e 
| *1 (9) = gantaytay + SST - tana ee 
: e 
| 
| ¥, (9) = tana . X, (9) (4. 26) 
re = ¥, (J) /sin [8 (J) + 6, (5) ] (4.28) 
X, + tana 
. sinT6 (3) + 6, (S77 - tana . cos [8 (5) + 567) 
(4.28.1) 


(b3) Knowing re, I,‘ (9) can 
be solved using equation (4.23). 

| - Calculate ve from eq. (4.22) 

| - Calculate E from eq. (4.2) 


(b4) Calculate second expression 


(b5) The iteration routine: 
The iteration search 
for correct Ss uses the Newton-Raphson routine. 

The function: 


Fer = 0 (4.29) 


ff ‘Ef 
is the expression from the formula (eq. 4.4) where r 
the geometric expression (eq. 2.28). 


for front radius (r¢) from eq. (4.4). 
The derivative is found with respect to e.: 


a(reg) d(r¢) 
4 —_ - 
P' = ate) ~ aay eee) 
d(rec) (r-.)3 es) 
ff arg: ff ‘ By ‘ cos (+20, + sin (Bo) ] (4.30.1) 
e 2E 
Where: 
2 2 
By = /n + 4B(E/r,) (4.30.2) 
2 2 
By = (2(E/r) -~ A]/B, (4.30. 3) 


The lower signs are used when the front radius (r-) is calculated 
from a known rear radius (x). (this case occurs at J > 1), while 


the upper signs are used when ry is calculated from ie (for Jwl). 


d(r,) -X, . tana. {cos [8, (J) + 6,(5)] + tana . sin[6)(J) + 6 (5) I} 


ate isin[@, (J) + 6,(3)] - tana . cos[@,(J) + @.(3))} 


(4.30.4) 
(c) Correct Solution 

After knowing the correct value 
for 8,09), the procedure used in (b) can be followed resulting 
in the correct values for: 
Ti(J), 8, X, (9), ¥,(9), re(J), I,'(), veld), E. Using 
equation (4.2.2) ve can be calculated. Equations (4.15), 
(4.16),(4.17) can be solved for n(J), I,*(3), B(J) as shown 


previously for the J=l1 case (equations 4.18, 4.15, 4.16). 


32 


4.3 Skew Rays 
4.3.1 Front Surface 
As in the homogeneous sace, the intercept 
point at the front surface (XO, YO, 20) can be found (equations 
Se19 ely 3519.2, .3+19%3)% 

Also, the direction cosines of the normal to the front 
surface can be calculated (equations 3.22,3.23). The direction 
cosines of the ray can be found as well using equation (3.26). 

4.3.2 Front Radius 
Knowing the intercept point at the front 


surface, the front radius can be calculated: 


re(3) = [(X0 = xX)? + (vo)? + (20) 714 (4.32) 


Knowing re(J), the initial index of refraction In, (r¢)] can 
be calculated using equation (4.3). 
4.3.3 Refraction at the Front Surface 
As for the homogeneous case (section 3.2.4), 
the refracted angle at the front surface can be calculated 
using equations (3.30), (3.31). 
4.3.4 Initial Direction Cosines for the Internal 
Refracted Ray 
The initial direction cosines of the internal 
refracted ray (K‘',L',M') are calculated the same way as for 
the homogeneous case (section 3.2.5, equation 3.32). The 


vector form of the initial direction of the ray is therefore: 


aA 


R=K'.e +L' .e,+M'.e, (4.32) 


x 
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4.3.5 Front Radius Direction Cosines 
The vector form of the front radius is 
defined as follows: 
> aA a a 
te = ley - Cy + Xgy + Cy + Fez * &z (4.33) 


Where the direction cosines are: 


ley = (XQ - XQ) /E (dS) 


Fey? S0Le es) (4.33,1) 


Yeo ZO/r _(J) 


4.3.6 Normal to Plane of the Curved Ray at (X0,YO,2Z0) 


The vector cross product of es and R is normal 


to the plane defined by the two vectors. 


ey ey ey 
NP, = Ip X R= © ex ley Teg (4.34) 
K' L' MoI 
NP, = NPey - Cy + NPey - Cy + NPeo . @o (4.34.1) 
Where: 
NP ey = ley e M‘ = Tez e L' 
NP ey = Xe - K' - Sex ° M' (4.34.2) 
NP eg s Tex e L' >= rey e K* 
oe np. = /wP,2 + nP_2 + NP,2 (4.34.3) 
f£ £X fY £Z peace 
The normalized direction cosines of the normal to the plane 
of the curved ray at (XO,Y0O,Z0) are therefore: 
NP ¢y/NP ¢ : NP ¢y/NP ¢ ; NP ¢7/NP (4.34.4) 
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Calculating the Angle Between the Ray and 
the Front Radius 


> 


Ze - R= [r,| . IRI . cos(v,) (4.35) 


Both the ray and the unit vector for the front surface radius 


are normalized, i.e. 


( 
= 


iFe| = Core + re + of ¥ 
(4.35.1) 
Fe (kt? + a'2 4 Mt) 2 


There fore: 


-1 
cos [rey - K' +24 Pa Te Keg ° M'] 


(4.35.2) 


£Y 


Calculating constant E 


Knowing Ver constant E can be calculated 


using equation (4.2). 


4.3.9 Intercept Point of the Curved Ray with Rear 
Surface (Point xXI,¥I,2I) 
4.3.9.1 Equations 
As previously, oe is the angle between 


the front radius and the rear radius. 


rex: Tex + Loy + Tey + Teg -Tez 7 cos [@, (J) ] (4.36) 


The direction cosines of the front radius were calculated in 
equation (4.33). Similiarly, the direction cosines of the 


rear radius are: 


= (XI - X/ r,(d) 


ry” ¥I/r (J) (4. 36.1) i 


ae 2I/r (J) 


Therefore, equation (4.36) can have the form: 


The rear radius is defined as: 


(xr = x)? + cyt)? + (21)? = Ce tay? (4.37) 


The rear surface was defined previously. Consequently: 
-(axt - b)? + (yt)? + (21)? = 0 (3.33) 


(XI - Xa) Fey + (YI) rey + (21) reo = r (J) - cos[6. (J) ] (4.36.2) 
a,b were defined by equations 3.33.1 and 3.33.2. 
Also, from geometry: 


(XI - Xo) = re (o) : cos[6) (J) ] (4.38) 


Another equation arrives from the connection between the two 
radii through eq. (4.4). 
4.3.9.2 Solution 
a. Assume 0, (9)- 
b. Calculate r(J) from the formula 
(equation 4.4). 
c. Calculate XI from equation 


(4.36.2): 


XI = r (J) P cos [@) (J) ] + Xo (4.36.3) 


d. Calculate radius from equations 


lb 


(4.37) and (4.36.3): 


tyr)? + (21) 2) 


2 
r.(J) [1 - cos [6,(9) 1] 


r.(d) ° sin[6)(J)] (4.37.1) 


a and b from equations (3.33.1) and (3.33.2) as in the 


homogeneous case. 
. 
f. Substituting into equation 
(3.33): 


XI = (R+b)/a = {r, (J) . sin[@,(J)] + b}/a (4.38) ; 


Radius = R 
e. Calculate parameters 
| where R is the radius defined in equation (4.37.1). 


g. Solving equations (4.36.3), 
(4.38) for 84 (J) one gets: 


po a OR 


’ P; 1,2 2 ’ ’ 
yt) = sin Ee + Vie" = cance (4.38) 
Where: 
A’ = (r( 17 (1 + a) 


2 (4.40) 
B' = [r(3) 1(b as Xq) 


2 2 
* [r.(9)] 


cl=(b-a.X)?-a 
h. Knowing @)(J), it is pos- 
sible to calculate new value for a, (9) and iterate through steps 


(a) - (g) until the difference in the values of 89 (J) is 


37 


| 


small enough. 


i. When an accurate value of 


84 (9) is received, XI is calculated from equation (4.36.3). 


YI and ZI are calculated using equations (4.36) and (4.37.1). 


Consequently: 


2 


_ BS (p")? = (at) (c*) 
YI = a WBN) AMC) 


ZI = v fx, (J) ‘ sin{@, (3) 1}? - {yz}? 


| 


(4.41) 


AO ee tel /legel- 


B" = {r (J) F cos[6, (J) ] - (XI - Xa) Fey) Fey/Teg 
(r,(9) 
{ 


c" = 


r 
. cos[@ (J) ] - (XI - Xgl? - jx,(o . sinfe,(a)1(? 
leg £Z 
(4.42) 


When lez = 0, the solution for YI and ZI is different: 


YI = {r (J) - cos[6,(J)] - (XI - Xo) Fex)/ley 


While ZI is calculated by equation (4.41). 


4.3.10 Direction Cosines of the Ray at the Rear 
Surface 


= a. Similarly to what was done at the front 


surface (equation 4.35): 


z.- R= [Z| . IR). cosy) (4.44) 


Both the ray and the rear radius are normalized; therefore: 


t t 
rex ‘ kK. + rey e L. + Teg M.. = cosy . (4.44.1) 


b. The ray remains within the same plane 
while transfering from the front surface to the rear surface. 
i.e.: 


‘ Cree a 
ex ° KL + NP ey ‘ LL! + NP co ‘ M,. = cos (1/2) 0 (4.45) 


c. The direction cosines of the ray are 


i] 2 i] 2 a 2 - 
(kK. )~ + (Ly )~ + (M »"~ =l (4.46) 
ad. Solution: 
(dl) NPp, # 0, NP, # 0, NPp, # 0 


SSS 
. vis") ~ 4(A") (C") 


SI 
a 

u 
vi 
to 

+ 


normalized, i.e. 


: 2 (NP py) (NP gg) 


(NP py) 


| (mM) = E- 2 (ue) oe 
| (K' 2) = [-NPgy - (L'L) = NPgg . (M'L) INP gy 
| Where: 

A = NPey - Tey ~ NPey - Try 

BS NP ges Epp = Wag + thy 

C = NPsy - cos(¥,) 

ee rae pe (4,48) 

£Xx 
2 

. Bi = (ae ) +1 


= (At) + (Bt) (2)? = ey 


— 2(A)(c)(B')_ _ (C') (Cc) 
(B) - ey 


= (B') (c\e ra 


In (4.47) the (+) sign is used when rey < 0 while the (-) 
sign is used when rey 2 0- 
(d2) In the above case, when B = 0 (look for 


definition of B in equation (4.48)) the solution for direction 


cosines of the ray is different. 


Cc 


A 
(AONE ae irn. at yn 
(OU) eg] 8 A) Et) 


(4.49) 


All the symbols used here were defined in equation (4.48) and 


the + sign is used as before. 


(d3) NP ey = 0, NP ey = 0, NP ey 70, rey ¥ 0 


(M') = 0 


page ee 
pag Be VB) — = (a (0%) 
OSS Se a 


(4.50) 
1 ' 
(L' ) = es [cos (y) - Ty: K | 
The (+) sign is used here. 
While in this case: 
r 2 
areas (24) 
ry 
r 
B" = cos (p) zx 2 
(ry) (4.51) 
cos(W_) 
co =[——=] -1 
ry 
(a4) In the above case, when rey = 0 the 
solution is: 
(MY) = 0 
} cos (4) 
el eae 
rx (4.52) 


2 
(BL!) = t/1 = (Rt) 
Again, the (+) sign is used here. 


(a5) NP = 0, NP cy #0, NP eo = 0, re2 #7 0 


£X 


The solution is similar to that in case (d3). 


! 
! 
ql 


(Lt) = 0 
B" . Jip")* = (at) (c*) 
ct) = Bye VOBT) | fat con) (4.53) 


= 
tt 


K'] 


1 
lap [cos (¥,) = rex ° r 


(A"), (B"), (C") are the same as defined in (4.51) except that 


xr is replaced everywhere by reg: The (+) sign is used here. 


ry 
(d6) In the above case, when ae 0 the 


solution is similar to that in case (d4): 


(4.54) 


a ey a 
tyl- (KT) 


® 
W 


Also here, the (+) sign is used. 


(47) NP, = 0, NPpy #0, NPgo # 0, ry = 0 


£X x 


NP 
eek > rer (4.55) 
(MS) wee (E'Y) 


’ - : - ' 2 = e 2 
(K'.) = evil (L' .) (M' .) 


. 
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The (+) sign is used here. 


(a8) As in (d7), but: 


As in (a7), buts: rey # 0 


_ (B") , vip")? = 4(a") (c") 
2(A") 


As in case (dl), the (+) sign is used when: < 0 while 


the (-) sign is used when: Try 


5. FLOW CHART OF PROGRAM LENS 
(Homogeneous Case) 


Note: Subroutines SHAPET, SKEWT compute the shape of lens 
and the skew rays as described in the previous sections (3,4). 
However, subroutines SHAPEC, SKEWC compute the shape and the 


skew rays of a different design of the lens. In this design, 


the fixed cone is the rear surface instead of the front | 
surface (compare figure 9 to figure 2.). The equations for 
this different design are not included in the current report. 
The equations are part of Capt. Carr's thesis [15] and will 


be described there. 
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Program LENS ~ Main 
Start 


Change 


YES 


B (J) 


NOSE 
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: 
/ I,R,T(1) 


i 


Ny No ON 


, TCL) © Tt) 
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ene wer Oe 


*) 


In program GRIN, which uses the same subroutine, here, 
data for XC, DIV is given. 
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"(3) 


I,(J) 


N 


Compute 


Compute 


N 
H 
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RI Um oe 5 


Sidedheeteeiacememneteatnass temeeabnal bansen deltaic athe 


Lae aa ET ye A ine cee em 
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Compute 
x, (1) 


Compute 
yo (1) 


Compute 
8 (1) 


Compute 


y, (9) 


y, (9) 


— 


Compute 


T(J) 


| Compute ! 


' 


x5 (3) py, (9) 


i Compute i 
| 
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ee | 
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Compute 


xy (J) r¥, (9) 


Compute 


I, (9) I," (g) 


Compute 


¥, (9) »x, (9) 


j 57 


Swe mere ee 


SSeS cee 


Shift: 


x, (J) , x, (9) 


Compute 

Statna 

Compute 
Y 


Compute 


Onp 


(*) In program GRIN, which uses the same subroutine, here, 
/ also xc is shifted. 58 


Ceni,Torbi, 
Boti 


Compute 


yg (12) 


61 


NzezeN-]1 


62 


Compute 


Cen, Torb 


Zo (12) = 0. 


ICHECK=60 


YE 
0 (12) <Bot 
NO 


ICHECK#20 


65 


66 


BA er ie A 


NO 


Return 


YES 


68 


=cos (0) cos (d) 
y 


=sin(d) cos (@) N_=-cos (a) 
2 2 


Compute 


xI, yI, 2I 


PRINT 


ICHECK=51 


tee are 


Compute 


{ b 
You (12) Zou I 


Compu te 
xDiapt (Iz) 


xDiapt (Iz) #10 


yDiapt (Iz) yDiapt (Iz) =10 


Flag (Iz)=0 


De fine 


¥q (12) ,24 (Tz) 


73 


ICHECK=0 
Nenety (Iz) =0. 


Flag(Iz)=1 


ICHECK=0 
Ntncty (Iz) =0. 


- Count + 1 


Count 


ICHECK=0 


tncty (Iz) =0. 


Flag(Iz)=3 


Flag (Iz) 


=Flag(Iz.) 


Count=Count 
+] 


) 
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ydiapt (Iz) 


sy diapt(Iz.) 


Opl (Iz) 
=0pl (Iz.) 
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Se Se 


Rays*Ia-1 
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Compute 
Ely(), Ells 


Coord=Coord 
+1 
El, Q =El. +EL 
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El OsE1,. 
-Ell 


Coor() =Coord 


Coord =Coord 
+1 


De fine 
Number ,Numb 


=e, .10 
P 


I 
angle 
Thick=T (3) 


Compute 
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Compute 
Sdrad (Iz) 
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Compute 
Fractn (Iz,) 


Roc(Iz,) 


=Roc (Iz, ) 
+0.005 
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PRINT (83) 


PRINT (84) 


PRINT /86) 
(87) 
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PRINT (89) 
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6. FLOW CHART OF PROGRAM GRIN 


Note: Program GRIN is identical to program LENS (Section 5) 
except for subroutines: 

SHAPET, SKEWT, SHAPEC, SKEWC. 
The flow chart for the first two is listed below. Subroutines 
SHAPEC, SKEWC are of a different design (fixed cone a, for 
the rear surface). They are part of Captain Carr's thesis [15], 


and they will be described there. 
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pea ee a EE ew oceeeiemiees ae 


eee Sette AEN EE AE RR rc sb ec tren a een er a ETE Oye 


wee ee pe er en er ry eS 


Kitter=1 
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t 


Compute 


= 2) 


moe 


Compute 
y,(1) Xo (1) 


+ ln aroha nom as 


PIER tr ee ee on 


Compute 


xy (J) r¥, (9) 


Compute 


re(d) 


IRAE IAN ENP 6 TTR NCRMORE Hr mn te 


Compute 


x0,y0,z0 


N =-sin (a) 


<> YES YES 
O 


N_=cos(¢) N =0 N =O. 
¥ cos (a) y f 
N_=sin(0) N=cos (a) N =-c0s (a) 


Change sign 


in Ny Ny 


o= 
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Compute 


vex!" ty'* tz 


Compute 
(E08) Wek 


8,3 (8,48) 


R=Radius= 


x -sin(8)) 


rf 


piff= 


o (NEW) 
89 (OLD) | 


NO 


Compute 
aa nr tad 


ICHECK=52 


wee 


PER PND 0 Sere rR or wetter --n- a ore ee nae 


Compute 


Return 


102 


Lt A a onerpene annie meme ee 


reg oy 


Compute 
¥ (6/8) ,Oxr 


re(x,.8) 
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Compute 


r,(J) -geom. 


DIRECT 
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7. COMPUTER PROGRAMS USER'S GUIDE 


Input Data 
F, a, U, I case 
I, R, T(1) 

Ni, No, N3 

ay Grid, Ellnum 

Xor Div 

Shape, Squray, Ellips 
Options 


I 1 Rear surface is a fixed cone 


case~ 
(Capt. Carr's thesis, see ref. [15]) 


I 2 Front surface is a fixed cone 


case 
(This report). 

Execution Commands 

7.3.1 Opening Commands 

LAXXXX (Write account no. instead of XXXX). 

Password 

GLOBAL TXTLIB FORTMOD2 MOD2EEH 

Compilation 

FORTGI LENS (or GRIN) 


Run on Terminal (worke‘)) 


FILEDEF 05 DISK DATA L 

FILEDEF 06 DISK OUT D (RECFM F BLOCK 132 PERM 
EXEC RUN LENS (or GRIN) 

XEDIT OUT D 


Hard Copy 
FILE 


PRINT OUT D 


Get a Plot 


LENSCOM PRINTER 
| 7.3.6 Change in Plot featine 
XEDIT LENSPLOT FORTRAN 
Transfer Files to Another User 


(a) 


The Transferer 


CP SPOOL PUNCH ConT TO xxxxp ‘?) 
PUNCH LENS FORTRaN ‘2? 
(PUNCH CRIN FORTRAN '2)) 


(puNcH pata 1 ‘2)) 


CP SPOOL PUNCH CLOSE NOCONT 
CP SPOOL PUNCH OFF 
The Receiver 
READCARDA *A* 
Notes 


(1) 


WORKG is the name of the file of 


commands in this section. 


(2) Account No. of the receiver. 


(3) Optional. 


8. COMPUTER PROGRAMS LIST OF SYMBOLS AND SUBROUTINES 


8.1 Subroutines 


DATTA Initial data 

SHAPET Computes shape of lens in a 

SHAPEC Computes shape of lens in case”! 

INIT Computes behavior of first ray. 

POINTS Computes coordinates of rear and front surfaces. 

FRONT Computes I, (J), I,'(J),¥,E at front surface. 

REAR Computes n(J) ,I5(9) ,I,' (J) at rear surface. 

CALTC Computes thickness at each point for SHAPEC. 

NOSE Computes extrapolation of last points in Lcagautt 
or coordinates shift in I ase ** 

ZOYO Computes Zo‘ (IZ) sy! (12) for skew rays calculation. 

RADI Computes intersection of the ray with the designed 
surface. 

SKEWT Skew rays calculation in E ekg e* 

SKEWC Skew rays calculation in L eea* 

INTENS Calculates transmitted intensity of the rays. 

IMAGE Calculates image plane spot diagram. 

SUMM Calculates spot diagram energy cGensity. 

PLOTT Save output for plot routine. 

INDEX Computes index of refraction in the lens (No). 

RADN 2 Computes radius to rear (or front) surface from 
known radius to front (or rear) surface and 6,. 

ITER Iteration routine to find correct 36 in SHAPET. 

DIRECT Computes direction cosines of the ray at the rear surface. 
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| 
| 


8.2 Program Symbols 
8.2.1 Data Symbols 


F Focal length. | 
a The fixed surface cone half angle. ! 
U The incident ray offset angle. 
I Number of iterations. 
R Maximum radius of fixed surface cone. } 
T(1) Thickness of lens at the edge, | 
Ny /No.N, Indices of refraction. | 
a, Inclination of object plane. 
Grid Size of grid increment. 
Ellnum Number of ellipses. 
Ell z - coordinate increment for ellipse plot. 
Xo Location of center symmetry in GRIN. 
Div The maximum change in the index of refraction 
(No) 

8.2.2 Geometrical Symbols 
Xx (5) py, (9) Coordinates on front surface. 
X(T) ,y2 (9) Coordinates on rear surface. 
(J) Radius to front surface. 
ry (d) Radius to rear surface - 
XO, YO, 20 Intercept of the ray with the front surface - 
XI,YI,2Z1I Intercept of the ray with the rear surface - 
Yom’ 21m Intersection with the image plane - 


8.2.3 Direction Cosines Symbols 


NUON iN, Direction cosines of the front surface normal - 


Y 
114 


Direction cosines of the rear surface normal. 


Direction of cosines of the ray. 

Direction cosines of the internal refracted 
ray. 

Direction cosines of the front radius. 
Direction cosines of the normal to plane of 
the ray at front surface. 

Direction cosines of the rear radius. 
Direction cosines of the rear surface external 
refracted ray. 


Angle Symbols 


The angle symbols are defined in ficqures 3 and 8. 
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9. PROGRAM LENS - LISTING 
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11. PRELIMINARY RESULTS 


with the RMS Spot size for a GRIN lens. 


0.0057336 0.0069021 
0.027424 0.027646 
0.06277 0.06288 


0.113431 0.11216 


0.18265 0.18287 
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Table 1. Comparison of RMS Spot size for a homogeneous lens 


0.005286 
0.027327 
0.06304 
0.113245 
0.18277 


Note: DIV is the fractional variation of refractive index. 
a, is the angle of incoming parallel rays relative to lens 

axis; a, has dimensions of radians. The center of symmetry 
Xo was at -(.07 measured from the apex of the outer surface. 


The minus sign indicates that X, is external to lens. 
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SUMMARY OF RESULTS 
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SUMMARY OF RESULTS 


FOCAL LENGTH = 4.0000000 ALPHA = 0.7853979 U = 0.0 


> R =» 1.0000090 T = 0.0500000 


ALFAP = 0.2000000 GRID = 0.10N0000 STATNA = ~O0.44L7L GAMMA = 0.7850314 
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12. SUMMARY OF THE REPORT 


Lens design and ray tracing to determine optical 
aberrations requires efficient computer codes. A clearly 
structured computer code with efficient algorithms has been 
developed to design a lens incorporating gradient refractive 
index, GRIN. One surface of the lens, either the inside or 
outside surface, is a circular cone. 

The computer code includes lens design and calculation 
of optical aberration for a lens using a homogeneous material, 
i.e., for n equal to a constant. The homogeneous case is 
included for several reasons. The simpler case with n equal 
to a constant provides insight to the more complex case of 
variable n. Extensive computer codes require a check on the 
validity of the output. The homogeneous case provides a check. 
When the factor B is equal to zero in equation [2.9], the 
results for the two computer programs with and without GRIN 
should be identical. Needless to say, the results are identical 
for the programs reported here. 

The case of a lens with GRIN and with a circular cone for 
the outside surface has been programmed. The companion case 
with the circular cone on the inside is being developed by 
Carr [15]. Carr's program can be readily incorporated into 
this computer code. 

The program has been exercised for the case of 5% 


variation in refractive index for both negative and pos.tive 
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values for B. The aim of the preliminary calculations 

was to verify operation of the computer program. More extensive 
calculations are required in order to understand the influence 
of various variables. As more experience is gained, the 
program should be combined with an optimization code to 

permit thorouch and efficient lens design. 

The selection of a 5% variation in refractive index 
was based on current state-of-the-art in the production of 
GRIN materials. A wider variation in refractive index may 
be needed. 

The program is restricted currently to values of Xo 
satisfying the constraint X, < 0. Additional work is needed to 
permit designs using positive values of Xo: When X, is 
positive, the angle © at some point along the ray interior to 
the lens has a critical value at which dr/d6 changes sign. 


The critical value causes computational difficulties. 
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